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00 '. Abstract 

' In this paper I show that it is possible to use Regge theory to constrain the 

I initial parton distribution functions of a global DGLAP fit. In this approach, 

OO . both quarks and gluons have the same high-energy behaviour which may also be 

I used to describe soft interactions. More precisely, I show that, if we parametrise 

■ the parton distributions with a triple-pole pomeron, i.e. like log^(l/x) at small 
I X, at — Qq and evolve these distribution with the DGLAP equation, we can 

■ reproduce F|, F^, F^/F^, F^'^ and xF^'^ for > 12.5 GeV^. In this case, 
I we obtain a new leading-order global QCD fit with a Regge-compatible initial 

■ condition. 

i| I shall also show that it is possible to use Regge theory to extend the parton 

• distribution functions to small . This leads to a description of the structure 

^ I functions over the whole range based on Regge theory at low and on 

^ '. QCD at large Q^. 

^ I Finally, I shall argue that, at large Q^, the parton distribution functions 
obtained from DGLAP evolution and containing an essential singularity at j = 1 
can be approximated by a triple-pole pomeron behaviour. 



X 



1 Introduction 



About thirty years ago, Dokshitzer, Gribov, Lipatov, Altarelli and Parisi have shown 
[1] that quantum chromodynamics predicts a breakdown of Bjorken scaling in Deep 
Inelastic Scattering (DIS). Once the parton distribution functions are fixed at one 
initial scale = Qqj the DGLAP equation gives their evolution to larger values of 
Q^. Although the initial equation has only included QCD contributions at leading 
order (LO) in Os, the next-to-leading order (NLO) corrections are now known as 
well as the NNLO corrections. There exists a rather large number of global fits (e.g. 
[2, 3, 4, 5, 6]) using the DGLAP equation to reproduce the DIS data. The basic 
idea is to fix the initial parton distributions, not predicted by perturbative QCD 
(pQCD), and to evolve it in order to reproduce the experimental measurements as 
well as possible. The success of this type of analysis is often considered as one of 
the most important prediction of pQCD. However, it appears that this approach 
presents some problems. Firstly, even if the strong rise of F2 observed by HERA at 
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small X is well reproduced by the DGLAP evolution, this may seem surprising since 
the evolution generates an unphysical essential singularity which should be replaced 
by the BFKL small-x behaviour [7] which is not observed in the data and is unstable 
against NLO corrections [8]. In addition, since the initial parton distributions are 
not predicted by pQCD, we have to parametrise them. This introduces a large 
number of free parameters in the models. 

In this paper, we shall study the possibility to use Regge theory [9, 10, 11] to 
constrain the initial parton distributions used in QCD global fits. To motivate this 
approach, one can, for example, consider the MRST2002 initial conditions: at small 
X, we have 

xq{x,Ql) = Ax-°-^^, 
xg{x,Ql) = Sx-°-27 + Cx°°^ 

These singularities in x do not correspond to any singularity present in hadronic cross 
sections [12, 13, 14, 15] and, conversely, cross-section singularities are not present 
in parton distributions. There should therefore exist a mechanism explaining how 
the residues of these singularities in partonic distributions vanish when goes to 
zero, and how the residues of the singularities observed in the total cross sections 
vanish for nonzero . Such a mechanism is unknown and seems forbidden in Regge 
theory, hence a description of both total cross sections and partonic distributions 
with the same singularity structure seems necessary. 

In a previous work [16], we have shown that, if one considers only the small-x 
and large-Q^ domain, on can use a triple-pole pomeron (squared logarithm of x) to 
reproduce the low-Q^ data and evolve it using the DGLAP equation to obtain the 
high-Q^ measurements. More precisely, we have used parametrisations of the form 

Alog^il/x) + Blog{l/x) + C + D 

for initial quark and gluon distributions. Since this Regge-constrained parametrisa- 
tion does not extend to a:; = 1, we have used the GRV parton distributions at large 
X {x > .TRcggc ^ 0.15). As a consequence, only two quark distributions were needed 
(a flavour singlet coupled to gluons and a flavour non-singlet) and only the proton 
structure function was fitted. 

In this paper, we shall extend the triple-pole parametrisation up to x = 1 in order 
to obtain a global QCD fit compatible with Regge theory at small Q^. At first sight, 
if we want to replace the GRV parton distributions at large x, we can, for example, 
introduce powers of (1 — x). However, if we do so, we must not only concentrate 
on the Regge domain but also on the large-x experimental measurements, including 
other experiments like j*d scattering and the neutrino data. We shall fit all structure 
functions over the whole x range at all scales greater than Qq. This will allow us to 
extract the parton distribution functions, which will be parametrised in such a way 
that they agree with Regge theory at small values of x. 
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In this study, wc shall also consider the extension of the Regge initial parametri- 
sation to low-Q^ values. Using standard techniques, we shall show that it is possible 
to continue the initial parton distributions to small Q^. 

Finally, we shall show that, at large Q^, the parton distributions obtained from 
DGLAP evolution, containing an essential singularity, can be approximated by a 
squared logarithm of 1/x within estimated DGLAP uncertainties. This confirms 
the results obtained in [17] where we have shown that the residues of the triple-pole 
pomeron can be extracted from DGLAP evolution. 

We should also point out that a similar type of approach has already been used by 
Donnachie and Landshoff [18], and by Csernai, Jenkovszky, Kontros, Lengyel Magas 
and Paccanoni [19]. However, our approach is different from their work in a number 
of points. Firstly, we keep the full DGLAP evolution equation, without taking 
only the residue at the leading singularity as done by Donnachie and Landshoff. In 
addition, the work by Csernai et al. only reproduces which means that they 
only need two quark distributions: a flavour singlet and a flavour non-singlet, only 
the flavour-singlet distribution being parametrised using Regge theory. Ther present 
work extends these ideas and those of [16, 17] to provide a standard set^ of structure 
functions, which reproduces the data for all values of x and Q^. 

2 Fitted and evolved quantities 

If we want to extend the paxametrisation introduced in [16] up to x = 1, we cannot 
only restrict ourselves to - In order to have a good determination of the valence 
quarks and of the sea asymmetry, we also need to include other structure functions, 
measured in the large-x region. In this global fit, we thus include the following 
quantities: 

• The proton structure function [21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 

33, 34, 35, 36, 37, 38, 39]: this is by far the most important type of experimental 
data. Moreover, it is nearly the only one to contribute to the fit in the small-x 
or in the high-Q^ region. 

• The dcutcron structure function F2 [37, 38, 39, 40]: as we shall see, these data 
allow the determination of the sea asymmetry. Many points are available in 
the large- and middle-a; regions, where the sea asymmetry is expected to be 
large. 

• The neutrino structure functions and F^^ [42, 43, 44]: these data, in 
which most of the points are at large values of x, are important to fix the 
strange quark and the valence quark distributions. Note that the data consid- 
ered here are averaged over neutrinos and anti-neutrinos. 

• The F2 /F2 measurements [41]: these data constrain the valence quark distri- 
butions and the sea asymmetry. 

C code for this standard set is available at littp://lepton. thco.phys.ulg.ac.be/ soyez. 
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Once wc know which experiments are fitted, we must find which quantities need 
to be evolved. Since the range under consideration in global fits extends up 
to 30000 GeV^, we need to consider 5 quark flavours: u, d, s, c and b. In order 
to use DGLAP evolution, it is easier to perform linear combinations of the quark 
distributions. In our case, we shall use 6 fiavour-non-singlet distributions 

xuv = x{u — u), 

xdy = x{d — d), 

Ts = xiu+-d+), 

Ts = x{u+ + d+ -2s+), ^ ' 

Tis = x{u^ + d+ + s+ - 3c+), 

= + d+ + s+ + c+ - 46+), 

where = q + q. Note that since the proton does not contain constituent strange, 
charm or bottom valence quarks, we have s = s, c = c and b = b. At leading order, 
the evolution of each of these distributions is given by the DGLAP equation 
with the splitting xPqq{x). In addition to the non-singlet distributions, we have the 
singlet quark distribution 

S = + d+ + s+ + c+ + 6+) 
which evolves coupled to the gluon distribution G = xg, with the full splitting matrix 

fxPgg{x) 2nfXPqg{x)\ 
\XPgq{x) XPgglx) J ' 

We shall assume that for < Am^, the quark q does not enter into the evolution 
equations. 

If we invert the relations (1) and express the quark densities q'^ in terms of the 
evolved quantities, we obtain 

XU+ = -^(12S + 3r24 + 5Ti5 + 10r8 + 30T3), 
60 

xd+ = ^(i2s + 3r24 + 5ri5 + ior8-3or3), 
XS+ = ^(i2E + 3r24 + 5ri5-2or8), 
XC+ = ^(4s + r24-5ri5), 
Xb+ = ^(E-r24). 

5 

Now, we can of course write the structure functions considered here in terms 
of the parton distributions or in terms of the flavour-singlet and flavour-non-singlet 
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distributions^. If, for the sake of clarity, we include other quantities like the neutron 
structure function, this gives 

= f (n+ + c+) + |(d+ + .+ + 6+) 

= (22S + 3T24 - 5ri5 + 5Ts + ISTg ) , 

= y(d+ + c+) + |(«+ + s+ + 6+) 

= (22S + 3r24 - 5Ti5 + 5T8 - ISTg), 
, _ Fj + F^ 

= ^(22s + 3r24-5ri5 + 5r8), 

and for the neutrino structure functions 





= 2x(d + s + 6 + n + c), 


pun 


= 2a;(it + s + 6 + (i + c), 




= 2x{u + c + d + s + h), 


pun 


= 2x{d + c + u + s + h), 


XF^/ 


= 2x((i + s + 6 — « — c). 


xF^ 


= 2x{u + s + h — d — c), 


xFl^ 


= 2x{u + c — d — s — b), 




= 2x{d + c — u — s — b), 



If we average over proton and neutron targets, we obtain the neutrino-nucleon struc- 
ture functions^ 

xF^^ = x{uv + dv + s+ ~ c+ + b+), 
xif^ = x{uv + dv - s+ + C+ -b+). 



^At leading order, the quark coefficient functions are proportional to 5(1 — x) and the gluon 
coefficient function vanishes. 

^Neutrino experiments are often performed with heavy nuclei which means that the averaged 
structure function is measured. 
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We may finally average over neutrinos and anti-neutrinos, which leads to 

F/^ = x{u+ +d+ + S+ + C+), 

= s, 

xf/^ = x{uv + dv). 

3 Initial parametrisation 

If we want to perform a DGLAP evolution, we need to fix the parton distribu- 
tion functions at an initial scale Qq. Following the same ideas as in [16], we shall 
parametrise each quark distribution as the sum of a triple-pole pomeron term and 
an a2//-reggeon term. In addition, each distribution will be multiplied by a power 
of {1 — x), to ensure that the parametrisation extended to x = 1 goes to when 
X ^ 1. This leads to the following parametrisation 

xq{x,Ql) = [A,log\l/x) + Sglog(l/x) +Cg + D,x^] (1 - xf\ 

with^ q = uv, dv, Ug, dg, Sg, Cg and g. Fortunately, we can restrict many of the 35 
parameters introduced here: 

• First of all, the charm (bottom) distribution will be set to zero for < 4m^ 
{Q^ < 4ml). shall therefore take Qq < 4m^ so that we can set c{x, Qq) = 
and b{x, Qq) = 0. In other words, we have Ti5{x, Q^) = T,{x, Q^) for < iml 
and T24{x,Q^) = S(x,Q2) for Q2 < 4^2 

• The pomeron does not distinguish between quarks and anti-quarks. This 
means that the valence distributions uy and dy do not contain a pomeron 
term. 

• The pomeron, having vacuum quantum numbers, is insensitive to quark flavour. 
Thus, the only parameter through which the quark flavour may influence the 
pomeron is its mass. In other words, the couplings Ag, Bg and Cg are func- 
tions of Q^ and only. Consequently, the pomeron contributions to the Ug 
and dg densities are the same. Assuming that the strange mass is very small 
compared to the virtualitics Q'^ under consideration, we shall also take the 
same pomeron contribution^ in Sg- 

■^u — -^d — Ag = A, 
Bu = Bd = Bg = B, 
Cu = Cd = Cg = c. 

''The sea distribution qa is simply . 

®If we insert an overall factor in Ss, the fit naturally sets it to 1. 
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We shall assume that the reggeon, being mainly constituted of quarks, does 
not contribute to the gluon distribution. The parameter Dg will thus^ be set 
to 0. 

We know from [20] that, at large x, the following behaviour is stable with 
respect to DGLAP evolution 

S ~ (l-xf, 
(1 - xf+^ 



G 



log (t^) ' 



The denominator log(l — x) in the gluon distribution does not have a good 
behaviour at small x so we have not included it^. Nevertheless, we shall impose 

bu = bd = bs = b, 
bg = b+l. 

• If we look at the large-x data, we can sec that if we only use Dx^{l ~ x)^ 
for the valence quarks, the resulting distribution is too wide, or has a peak at 
too small a value of x. In order to solve that problem, we have multiplied the 
valence-quark distributions by a factor (1 + ^qx). 

• Finally, we still need to impose sumrules. Quark-number conservation can be 
used to fix the valence-quark normalisation factors. If we write 



k = h 



= , ( 1 + „ , r , J • (2) 



we find 

Tjbq + l)r(r?) / 7gr? 
r(r/ + 6g + l) V ri + bg + 1^ 

The momentum sumrule is used to fix the constant term Cg in the gluon 
distribution. Although all the functions involved are analytically integrable, 
the resulting expression for Cg is quite complicated and we give it in appendix. 



*If we do not impose Dg — 0, the parameter stays small in the fit. 

^Onc solution is to multiply the gluon distribution by an overall factor x. This makes no change 
at large x and ensures a good behaviour at small x because, when a; — > 0, 

1. 



Numerically, including this factor in the gluon distribution only makes a small correction. 
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Taking all these considerations into account, we obtain the following parametri- 
sation for the initial distributions 



XUy 


= l^x^{l+j,x){l-xf-, 


xdy 


= ^x\l+^dx){l-xf\ 

d 


XUg 


= [A\og^{l/x) + Slog(l/x) + C + Dux"^] (1 - xf 


xds 


= [A\og^{l/x) + Slog(l/x) + C + Ddx"^] (1 - xf, 


XSg 


= [A\og^{l/x) + Slog(l/x) + C + Dsx"^] (1 - xf, 


XCg 


= 0, 


xbs 


= 0, 


xg 


= [Ag\og\l/x) + B,log(l/x) + C;] (1 - xf+\ 



(3) 



where the parameters marked with an asterisk are constrained by sumrules. 



4 Fitted experiments 

As said previously, we have fitted F|, F^, F^^ , xF^^ and F^/Fl^. We shall now 
detail which experiments are included in the fit for all these quantities. 

For the proton structure function, we have fitted the experiments from^ HI 
[21, 22, 23, 24, 25, 26, 27], ZEUS [28, 29, 30, 31, 32, 33, 34, 35], BCDMS [36], E665 
[37], NMC [38] and SLAC [39]. For the deuteron structure function measurements, 
we have included data from BCDMS [40] E665 [37] and NMC [38]. We have also 
taken into account the measurements of F2 /i^f from NMC [41] . Finally, the neutrino 
data used here come from CCFR [42, 43, 44]. 

Among all these experimental papers, some give, besides the statistical and the 
systematic errors, an additional normalisation uncertainty. For each of these subsets 
of the data, we have allowed an overall normalisation factor. Let Ri be the normal- 
isation uncertainty for the subset and pi the effective normalisation factor. We 
may easily minimise the with respect to this parameter by requiring 

dpi dpi ^ e^- 

where j runs overs the data in the subset i, dj^ ej and tj are respectively the jth 
data, its uncertainty and the associated theoretical prediction. We easily find 



i ^2 



Pi 



'3 7 



J2i 7^ 



*The dataset is coming from the DURHAM database (http://durpdg.dur.ac.uk) to which we 
have added the 2000 and 2001 data from HERA [26, 27, 35] as well as the reanalysed CCFR 2001 
data [44]. 
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Parameter 


Value 


Error 


A 


0.00876 


0.00043 


B 


0.0197 


0.0035 


C 


0.000 


0.017 


An 

y 


0.258 


0.032 


Bg 


-0.62 


0.25 


Du 


0.378 


0.030 


Dd 


0.480 


0.030 


Ds 


0.000 


0.013 


V 


0.392 


0.019 


lu 


7.46 


0.91 


Id 


9.1 


1.6 


K 


3.625 


0.016 


bd 


5.261 


0.086 


b 


6.67 


0.27 


Nu 


2.015 




Nd 


1.723 




Cy 


3.158 





Table 1: Values of the fitted parameters in the parton distributions. The last three 
parameters are not fitted but are obtained from sum-rules. 

Finally, we shall require that pi does not lead to a normalisation bigger than the 
uncertainty Ri. This means that we shall constrain pi to verify 

1- Ri< Pi<l + Ri. 

Before going to the result, one must point out that we have used here the latest 
CCFR data^ from 2001 [44]. These data from U.K. Yang's thesis are used by adding 
the errors in quadrature and, in order to solve a discrepancy with the other data, 
we have also allowed an overall normalisation factor of at most 3%. 



They consist into a reanalysis of the 1997 data. 



Experiment information 


This fit 


CTEQ6 LO 


CTEQ6 NLO 


Quant. 


Colab. 


Reference 


Nb Pts 


Pi (%) 


x' 


X^/nop 


norm. 


-1 norm. 


norm. 


-1 norm. 




BCDMS 


PLB223(1989)485 


167 




154.607 


0.926 


5.303 


5.303 


2.652 


2.652 




EG65 


PRD5i(199G)3(J0G 


30 


1.80 


i0.3G8 


1.316 


1.177 


1.233 


1.251 


1.383 




HI 


EPJC19(2001)269 


126 


-1.50 


129.673 


1.029 


1.516 


1.626 


1.077 


1.122 






EPJC21(2001)33 


86 




75.774 


0.881 


0.942 


0.942 


1.008 


1.008 






EPJC13(2000)609 


130 


-1.50 


117.682 


0.905 


1.612 


1.962 


0.882 


1.032 






NPB47()(199G)3 


15G 




104. 20G 


0.GG8 


0.835 


0.835 


0.G58 


0.G58 






NPB439(1995)471 


90 


-4.50 


49.499 


0.550 


0.597 


0.901 


0.574 


0.737 






NPB407fl993)515 


21 


-8.00 


6.233 


0.297 


0.289 


0.466 


0.287 


0.401 




NMC 


NPB483(1997)3 


79 


2.10 


101.927 


1.290 


1.728 


1.260 


1.138 


1.186 




SLiVC 


PLB282(1992)475 


52 




97.8G1 


1.882 


2.123 


2.123 


1.355 


1.355 




ZEUS 


EPJC21(2001)443 


214 




207.294 


0.969 


2.454 


2.454 


0.875 


0.875 






EPJC7(1999)609 


12 




11.297 


0.941 


0.744 


0.744 


1.259 


1.259 






ZPC72(1996)399 


172 




238.882 


1.389 


1.299 


1.299 


1.429 


1.429 






ZPC65(1995)379 


56 


2.00 


27.477 


0.491 


0.495 


0.415 


0.453 


0.470 






ZPC69(1995)607 


9 


-1.54 


11.493 


1.277 


1.201 


1.270 


1.309 


1.289 






PLB316(1993)412 


17 


6.94 


6.048 


0.356 


0.370 


0.372 


0.344 


0.474 




Total 




1417 




1380.321 


0.974 


1.864 


1.900 


1.150 


1.187 




BCDMS 


PLB237(1989)592 


154 




127.941 


0.831 


1.546 


1.546 


0.903 


0.903 




E665 


PRD54(1996)3006 


30 


- 


33.563 


1.119 


0.913 


0.913 


1.132 


1.132 




NMC 


NPB483(1997)3 


79 


1.00 


90.330 


1.143 


1.438 


1.131 


0.969 


1.071 




SLAC 


SLAC-357(1990) 


50 




98.376 


1.968 


2.515 


2.515 


1.278 


1.278 




Total 




313 




350.210 


1.119 


1.613 


1.536 


1.002 


1.027 


Fr 


CCFR 


UK. Yang's thesis 


65 


3.00 


165.512 


2.546 


3.118 


4.570 


3.523 


6.135 


xF^^'^ 


CCFR 


PRL79(1997)1213 


76 




42.066 


0.554 


0.658 


0.658 


1.252 


1.252 


F^/F^ 


NMC 


NPB371(1995)3 


91 




116.720 


1.283 


1.315 


1.315 


1.285 


1.285 


Total 






1962 




2054.830 


1.047 


1.794 


1.855 


1.215 


1.333 



Table 2: Fit results detailed experiment by experiment. For comparison we have added the predictions for CTEQ6 at leading 
and next-to- leading order (the NLO predictions are taken in the DIS scheme). In the comparison with CTEQ, the results are 
given with and without taking into account our normalisation factors. 
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5 Results of the DGLAP global fit 

We have adjusted the 14 parameters A, B, C, Ag, Bg, D^, D^, Dg, b, bu, ba, 7u, 7d 
and r) to the experimental data in the region 

q2 > = 6.76GeV2, 

> 12.5 GeV^. 

The second boundary is used to cut the region where higher-twists effects are ex- 
pected to be large and we have adopted the same limit on W'^ as MRST. The values 
of the fitted parameters are presented in Table 1 and the result, detailed experiment 
by experiment, is given in Table 2. In addition, the curves resulting from our fit are 
presented for each experiment in Figures 3 to 15. 

We can see from the parameter table that both the large-x exponents and the 
reggeon intercept have acceptable values. 

In order to evaluate the quality of our fit, we have also shown in Table 2 the 
CTEQ6 results at LO and at NLO (in the DIS scheme^'^), with and without taking 
the normalisation factors into account^^. We see that the CCFR 2001 neutrino data 
probably need to be renormalised up and are still poorly reproduced. We can also 
see that, apart from the SLAG data, we obtain a very good description. This means 
that it would be a good idea to add a renormalisation factor of a few percents to 
the SLAG F| and data. 

The correlation matrix for the parameters is presented in Table 5. 

In Figure 1, we have shown some typical distributions and their evolution. 
The xuy and xdy valence quarks distributions both present a peak around x ~ 
0.1 — 0.2 and are, roughly speaking, within a factor 2. The sea asymmetry d — u can 
be written in the following form: 

, J xuv — xdv — T3 
x{d -u) = 

= (L»d-L»„)x''(l-x)^ 
This distribution has a maximum for 

0.1 for xuvi 
0.07 for xdv, 
0.056 for x{d — u). 

The evolution in of these three distribution shows the same behaviour: the peak 
is moved to smaller values of x and tamed while its width grows. We have also 
shown in Figure 1 the gluon distribution which grows quickly with . 

^°The DIS scheme is the renormahsation scheme where, at any order, the quark coefficient function 
is (5(1 — x) and the gluon coefficient function vanishes. 

^^The CTEQ6 results are obtained by using the last CTEQ parton distributions to predict struc- 
ture functions without any refit. Therefore, the LO and NLO results are just given for comparison. 
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The parton densities at various scales are plotted in Figure 2. First of all, when 
= Qo = 4m^, wc have no charm or bottom quark and both quark and gluon 
distributions are described by Regge theory, more precisely by a triple-pole and a 
reggeon contribution. At higher virtualities, charm quarks are non-vanishing and, 
for > Ami, ^^^^ have b quarks. For > Qq, the parton distributions have 
an essential singularity at j = 1. 

Finally, we can estimate the uncertainty on the initial distributions in the fol- 
lowing way: for the sea quarks or for the gluon, we have {D = for the gluon 
distribution) 

xq= [A\og^{l/x) + Blog{l/x) + C + Dx"^] (1 - x)^ 

If we assume that the uncertainties on the parameters are uncorrected, we obtain 
easily 

{Sxqf = {\og^{l/x)SA'^ + \og'^{l/x)dB'^ + SC'^ + [SD'^ + log\l /x)D'^Sr]'^] x"^" 
+ [Alog^il/x) +B\og{l/x) + C + Dx"^] log2(l - x)6b'^} (1 - xfK 
For the valence quarks, the initial distribution has the form 

xqv = Kx'^il - x)''{l + -fx) 

with K fixed by quark number conservation, and wc find that the uncertainty is 

{5xqf = K^x'^'^il - xf" { [\og^{l/x)8rf + log2(l - x)6b^] (1 + -fxf + x^d-f'^] . 

The resulting uncertainties on the initial distributions are shown in Figure 16, 
where we have also plotted the uncertainties obtained by taking into account corre- 
lations between the parameters (see Table 5). We see that this "traditional" way of 
estimating errors leads to much smaller uncertainties than the joint consideration of 
forward and backward evolution obtained in [17]. 

6 Regge theory at low 
6.1 Motivation 

If DGLAP evolution gives the behaviour of the parton distributions at large Q^, 
we expect soft physics to be described by Regge theory. In other words, Regge 
theory should not only be able to describe the initial DGLAP condition but also the 
structure functions for < < <3o- -'■'^ this section, we shall therefore try to extend 
the parton distribution functions at low values of Q^. Note that in this region, we 
cannot use the DGLAP equation anymore. In addition, if we want to use Regge 
theory, we must still restrict ourselves to the high-energy domain. Hence we keep 
the constraint 

W'^ > 12.5 GeV^ 

but allow to be in the region 

< Q2 < 6.76 GeV^. 
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Experiment information 


This fit 


Quant. 


Colab. 


Reference 


Nb Pts 


Norm. 




/nop 




E665 


PRD54(1996)3006 


61 


1.80 


55.265 


0.906 




HI 


NPB439(1995)471 


3 


-4.50 


0.834 


0.278 






NPB470(1996)3 


37 


_ 


15.518 


0.419 






NPB497(1996)3 


44 


-3.00 


35.660 


0.810 






EPJC21(2001)33 


47 


- 


61.721 


1.313 




NMC 


NPB483(1997)3 


67 


2.10 


46.473 


0.694 




SLAC 


PLB282(1992)475 


94 




102.417 


1.090 




ZEUS 


ZPC69(1995)607 


14 


1.54 


29.452 


2.104 






ZPC72(1996)399 


16 




11.339 


0.709 






PLB407(1997)432 


34 




10.753 


0.316 






EPJC7(1999)609 


32 




34.541 


1.079 






EPJC12(2000)35 


70 




86.120 


1.230 






EPJ21 (2001)443 


28 




48.034 


1.716 




Total 




547 




538.127 


0.984 




E665 


PRD54(1996)3006 


61 




76.977 


1.262 




NMC 


NPB483(1997)3 


67 


1.00 


40.064 


0.598 




SLAC 


SLAC-357(1990) 


98 




86.511 


0.883 






PRD49(1994)5641 


1 


0.89 


0.000 


0.000 




Tot al 




227 




203.552 


0.897 




CCFR 


UK. Yang's thesis 


19 


3.00 


63.159 


3.325 




CCFR 


PRL79(1997)1213 


35 




39.201 


1.120 




NMC 


NPB371(1995)3 


120 




101.445 


0.845 


Total 






948 




945.485 


0.997 



Table 3: Result of the small-Q^ fit detailed experiment by experiment. 
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Parameter 


Value 


Error 


aB 


15.05 


1.80 


Qa 


6.37 


2.72 


^ U 


1.885 


0.522 


Ql 


10.00 


8.70 




0.437 


1.07 




5.19 


3.11 




3.64 


1.87 




5.10 


2.38 


Ql 


87.9 


15.8 


£A 


1.002 


0.328 


£B 


0.581 


0.110 


£C 


18.19 


4.37 


^■Du 


0.340 


0.157 


^T>d 


1.618 


0.615 


£b 


1.916 


0.453 


^bu 


2.558 


0.593 


^bd 


10.00 


6.94 



Table 4: Value of the parameters with their errors for the low-Q^ fit. The scales are 
given in GeV^. 

6.2 Small-Q^ parametrisation 

If we want to use Regge theory in the small-Q^ region, we need to parametrise the 
parton distribution functions. We shall use the same expressions as in (3) with an 
additional Q^ dependence. However, if we want to consider the extension down to 
= 0, we know that we should use the Regge variable ^ = ^ instead of x. This 
means that we shall use the following distributions: 



xuviy, Q^) 
xdv{i^, Q^) 

XUs{v, Q^) 
xds{v, Q"^) 



N' 



-{2v)-^ 



2u 



Q_ 
2v 



2v ) 



(4) 



{AiQ^) [log{2,y)-B{Q^f + C{Q^)+Vu{Q'){2u)-^] (l 
{AiQ^) [log(2i.) - ^(Q2)] 2 + c(g2) + p,(g2)(2i.)-'?} 
{AiQ^) [log(2i.) - ^(Q2)] 2 + c(g2) + p,(q2)(2^.)-'?} 



Q2.biQ^) 

2v ) 
2v J 

2\ 



1 - ^ 



Q 
2u 
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where, once again, and are constrained by quark number conservation. We 
shall require that the parameters in these distributions match the initial distribution 
taken for DGLAP evolution at 6.76 GeV^. Using parametrisations of the form^^ 




li{Q^)=l^{Ql) 
P.s =0 



for (j) = A,C,Vu,Vd,b,bu,bd 



for i = u,d, 



and constraining the parameters a_A, a^, ac, ax>^, aj)^, ax>^, ab, a^^ and a;,^ with the 
DGLAP initial condition at = Qg ~ ^-^^ GeV^, we are left with 17 parameters: 
OB, Q\, Qg, Q^, Ql,^, Q\)^, Ql, Ql^, Ql^, £b, £c, £r>„, £r>d' ^b, ^ba- The 
expressions obtained once the constrained have been imposed are the following: 

Qo + Qa^ ^ 



buiO") 
bdiQ^) 



A * 

Ql 



as 



Q2 x=« 



Ql 



c 



4^/ Ql \Q'' + Qc 



Ql + Qb 
Qo "I" Qc \ 



"Ql 
Q' 



{Qir 



D^^^iQir 



Qo + Qt>^ 
Q' + Qk. 
Ql + Qvd 



= 



b 



0! 

'Ql 
'""Ql* 
''Ql' 
Id 



Ql + Ql''" 



q^+q: 



Ql + Ql 



Q^ + Q^ 



' Ql + Ql 



£B-\ 



B 



^^Since Vg already vanishes at = Qo, we have set it to zero in the whole small-Q^ region. 
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6.3 Dataset and systematic errors 

In the small-Q2 region {W^ > 12.5 GeV^, < 6.76 GeV^), we shall fit the same 
quantities as previously and the data coming from the same collaborations: 

• Ff: HI [22, 23, 24, 27], ZEUS [30, 31, 32, 33, 34, 35], NMC [38], E665 [37], 
SLAG [39], 

• F^: BGDMS, NMG [38], E665 [37], SLAG [45], 

• F^/F^: NMC [41], 

• F^^ and xF^^: GGFR [44]. 

Concerning the treatment of the systematic errors, we have used the same correc- 
tion factors as the ones obtained in the DGLAP global fit for the papers containing 
data in both the small- and the large-Q^ region and leave this factor free for the 
papers containing only data at small Q^. 

6.4 Results 

The parametrisations described above have been fitted to the 948 data in the small- 
region using MINUIT. The results of this fit, experiments by experiments, to- 
gether with the parameter values is presented in Tables 3 and 4. We see that, apart 
from the ZEUS 1995 data and the GGFR F2 data, we obtain a good description of 
the structure functions in the low-Q^ region. The per data point is quite good, 
considering that we have applied Regge theory to quite a large region as compared 
to usual approaches [13, 46, 47, 15, 48, 49]. The poor description of the ZEUS 1995 
data may come from the fact that the systematic uncertainties have been fixed in 
the DGLAP fit. 

The results for the structure functions in the low-Q^ region are shown in Figures 
3 to 15 together with the large-Q^ results. In the small-Q^ region, the curves are 
only drawn in the fitted region {W^ > 12.5 GeV^). We see that the experimental 
measurements are well reproduced. 

Finally, the curves in Figure 2 show the parton distribution functions obtained 
at small . It is interesting to notice that valence quarks are large at small and 
12.5 GeV^ which, in this case, correspond to small values of x. For example 
at Q2 = 0.1 GeV^ and x = 0.008, we have xu 0.043, xdsi 0.053, xuy «i 0.057 and 
xdv 0.019. 

7 DGLAP vs. Regge at high 
7.1 Motivation 

As we have shown in [17], we can consider that Regge theory also applies at large 
Q^. In these conditions, since Q^-dependent singularities are forbidden in Regge 
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theory, we expect a triple-pole behaviour at all values of Q^. The unphysical essen- 
tial singularity generated by DGLAP evolution should therefore be considered as a 
numerical approximation to a triple-pole pomeron at small x. 

Given these considerations, we have shown [17] that, using both forward and 
backward evolution, it is possible to describe the small-x experimental data with 
parton distribution functions of the form 

A{Q^) log2(l/x) + B{Q^) log(l/x) + C(Q2) + D{Q^)x'' 

where the Q^-dependent couplings are extracted from the DGLAP evolution equa- 
tion. 

7.2 Parametrisation and uncertainties 

In this QCD global fit, we would like to test if it is still possible to consider the 
result of the evolution as an approximation to a triple-pole. To achieve this task, 
we shall fit the parton distribution functions at each value of with the following 
form 



xuv 


= A^'?(l + 7,x)(l-x)''^ 






xdy 








XUs 


= [A\og^{l/x) + B\og{l/x)+C + Dux'^\ 


1(1 




xds 


= [A\og^{l/x) + B\og{l/x) + C + Ddx'T 


1(1- 


-x)^ 


XSg 


= Ns[A\og^{l/x)+B\og{l/x) + C + Ds 


:X'>] 


(l-x)^ 


XCs 


= Nc[A\og^{l/x)+B\og{l/x) + C + Dc 


:X''] 


(l-x)^ 


xbs 


= Nb [A log2(l/x) + B log(l/x) + C + Dbxi] 


(l-x)", 


xg 


= [Ag log2(l/x) + Bg log(l/x) + C;] (1 - 




+1 



Before performing this fit, we shall determine an uncertainty on the parton dis- 
tribution. Since we want to show that the DGLAP evolution generates an essential 
singularity which mimics a triple-pole behaviour, we should estimate the error in- 
troduced by the evolution. Since we have used LO DGLAP evolution, we estimate 
that the errors are of the order of the NLO corrections: 

qNLoix, Q"^) (l± as{Q'^))qLo{x, Q^). 

If we require that the intimal parton distribution at = Qq remains fixed, this 
leads to 

norm/ /-)2\ ^ ^ ^ <^siQ ) ^ / 

1 ± as{Qo) 

or, keeping only the leading term in the strong coupling constant, 

Aq{x, Q2) = \as{Ql) - as{Q^)\ qLo{x, Q^). 
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In addition, we shall take into account the fact that, at small Q^, there may also 
be higher- twist corrections. If we assume^^ these are at 5% at = Qq, we shall 
finally consider 



[l-x)Q^ 



qLoix,Q^). 



Finally, additional powers of 1 — a; are expected to describe the large-a; behaviour 
of the parton distributions at large Q^, hence we shall only consider the region 
10~^ < X < 0.1 for the sea quarks and the gluons. For the case of the valence 
quarks, we hope that the factor (1 + 7a;) is sufficient to reproduce the distribution 
for 10"^ < X < 1. 



7.3 Results 

To perform the fit, we have taken the parton distributions with their estimated 
uncertainties in 80 points regularly spaced in log(a;). Since the parametrisations for 
the valence quarks and the sea quarks have disjoint parameters, we have performed 
two different fits at each Q^. The result of these are presented at Figure 8 for 
Q2 = 100 GeV^ and = lOOOO GeV^. In addition, the predictions for the triple- 
pole residues for Fj* are shown at figure 8 together with the per point of the 
fit. 

We clearly see that the parametrisation works very well for valence quarks at all 
values of Q^, as well as for the Ug and dg distributions. For heavy quarks and gluons, 
although the jnop remains less than 1, there are some discrepancies between the 
DGLAP essential singularity and the triple-pole fit at small x and large Q^. However, 
these differences are only present for > 3 TeV, which means that it is impossible 
to distinguish between the two approaches in present experimental measurements. 
This high-energy region should be reached at the LHC and should provide very 
useful information to distinguish between the different models. 

Moreover, we have also tried to estimate the uncertainties on the parton densities 
in another way. We have considered the uncertainties obtained from the DGLAP fit 
at Q = Qq* 

Ql) - Aq{x, Ql) < q{x, Ql) < q{x, QI) + Ag(x, QI). 

We can then evolve g(x, Qg) ^ A.q{x, Qq) and using these quantities to obtain the 
uncertainties at all values of Q^. If we do so, we obtain very similar conclusions. 



8 Conclusions and perspectives 

We have seen that we can use Regge theory to constrain the initial parton densities 
at = Qq and obtain the distributions at higher virtualities with the DGLAP 

^^This value reproduces errors comparable with those obtained from our uncertainties estimation. 

In addition, the higher-twist term is relevant for middle-range values of while the subleading 
corrections are important at large . 
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evolution equation. In this approach, Regge theory is used to describe the low-Q^ 
data and QCD appUes at large Q^. In such a way, the complex-j-plane singularities 
are common to parton distribution functions in the initial condition and to soft 
amplitudes which provides a unified description at high-energy in the soft region. 

We have also shown in this paper that it is possible to define the parton dis- 
tributions in the low-Q^ region and to parametrise them using Regge theory. This 
parametrisation is useful to describe the DIS structure functions but should be used 
with care. Actually, since factorisation is not proven at small Q^, we cannot ensure 
that the parton distributions can be extended to = 0. Using our parametrisation 
to describe processes such as jet production may be incorrect. 

Considering the low-Q^ parametrisation together with the Global QCD fit, we 
have a combined description of the hadronic structure functions over the whole 

range. This model, consistent with DGLAP evolution and with Regge theory, 
reproduces the experimental measurements with a very good x^- 

In addition, we extended the approach of [16] to x = 1 using only forward 
evolution. We have not applied the techniques developed in [17] and extracted 
the behaviour of the fitted parameters by combining forward and backward 
evolution. The reason is that, even with a few parameters, there often exist multiple 
minima and it is quite hard to obtain a continuous result for all parameters. This 
situation is expected to be even worse with the parametrisation used here due to 
the larger number of parameters. Hence, in order to test the compatibility between 
the DGLAP-evolved parton distributions and a triple-pole parametrisation, we have 
shown that the parton densities can be approximated by a log^(l/a:) behaviour at 
small X and large Q^. This approximation works very well up to ^/s 3 TeV, 
and at which point it deviates from DGLAP for the heavy quarks and the gluons. 
This means that wc expect high-energy corrections to be important in this domain 
and that the LHC should provide very useful information to distinguish between 
the different high-energy models. In this high-energy region, one should expect 
contributions from the BFKL equation as well as unitarity and saturation effects. 

Finally, a NLO analysis will be performed in the near future. This gives a much 
more reliable description of the data, allows a more complete comparison with other 
parametrisations and gives a description of the Fc and structure functions. 
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Figure 1: Typical momentum distributions inside the proton at various Q^: 
u valence quarks, (b) d valence quarks, (c) sea asymmetry d — u and (d) g] 
distribution. 
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Figure 2: Quark distributions inside the proton at various Q^: (a) = 0.1 GeV^, 
(b) Q2 = 1 GeV2, (c) Q"^ = Ql = Ami, (d) = 4^^, (e) = 2OOO GeV^ and 
ii\ = 30000 GeV^. See 6.2 for details concerninH: the oarton distributions at 
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Figure 3: DGLAP evolution results for BCDMS i?f data {i = for the upper curve 
and is increased by 1 from one curve to the next one). 
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Q^CGeV^) 



Figure 4: DGLAP evolution results for E665 data, z = for the upper curve 
and is increased by 1 from one curve to the next one. 
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Figure 5: DGLAP evolution results for NMC data (the SLAC data appearing m 
the NMC bins have been added to the plot), i = for the upper curve and is 
increased by 1 from one curve to the next one. 




Figure 6: DGLAP evolution results for HERA data {x < 0.001). i = for the 
upper curve and is increased by 1 from one curve to the next one. 




Figure 7: DGLAP evolution results for HERA i?f data (0.001 <x< 0.005). i = 
for the upper curve and is increased by 1 from one curve to the next one. 
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Figure 8: DGLAP evolution results for HERA i?f data (0.005 <x < 0.04). i = 
for the upper curve and is increased by 1 from one curve to the next one. 
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Figure 9: DGLAP evolution results for HERA i?f data (0.04 < x). i = for the 
upper curve and is increased by 1 from one curve to the next one. 
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Figure 10: DGLAP evolution results for BCDMS data. ? = for the upper curve 
and is increased by 1 from one curve to the next one. 
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Figure 11: DGLAP evolution results for E665 data, z = for the upper curve 
and is increased by 1 from one curve to the next one. 
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Figure 12: DGLAP evolution results for NMC data (the SLAC data appearing 
in the NMC bins have been added to the plot), i = for the upper curve and is 
increased by 1 from one curve to the next one. 
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Figure 13: DGLAP evolution results for CCFR F^^ data, i = for the upper curve 
and is increased by 1 from one curve to the next one. 
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Figure 14: DGLAP evolution results for CCFR xFg data, z = for the upper 
curve and is increased by 1 from one curve to the next one. 
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Figure 15: DGLAP evolution results for NMC F^/F^ data, i = for the upper 
curve and is increased by 1 from one curve to the next one. 




Figure 16: Initial distributions with their uncertainties: the dark region represents 
the correlated uncertainties while the light one is obtained without taking into ac- 
count the correlations between the parameters. 




Figure 17: Triple-pole pomeron fit to the parton distribution functions obtained 
from DGLAP evolution. The first column shows distributions at = 100 GeV^, 
the second corresponds to = 1000 GeV^ and the third to = 10000 GeV^. 
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Figure 18: Triple-pole form factors for at large presented together with the 
result of the fit of a triple pole to the large-Q^ parton densities/. 
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-0.494 


0.001 


0.521 


0.040 


0.005 


0.380 


-0.443 


0.900 


1.000 


0.081 


0.648 


0.552 


-0.893 


bu 


0.97469 


0.018 


-0.081 


0.000 


0.267 


0.397 


0.004 


0.253 


-0.317 


0.253 


0.081 


1.000 


-0.232 


0.537 


-0.073 


bd 


0.97329 


-0.022 


-0.042 


0.000 


0.085 


-0.243 


0.000 


0.199 


-0.228 


0.318 


0.648 


-0.232 


1.000 


0.174 


-0.332 


b 


0.99685 


0.442 


-0.589 


0.002 


0.876 


0.762 


0.008 


0.356 


-0.464 


0.518 


0.552 


0.537 


0.174 


1.000 


-0.389 


V 


0.99965 


-0.368 


0.562 


-0.001 


-0.460 


-0.009 


-0.004 


-0.291 


0.332 


-0.980 


-0.893 


-0.073 


-0.332 


-0.389 


1.000 



00 

o 



Table 5: Fitted parameters correlation coefficients. 



CO 

00 



A MOMENTUM SUM RULE AND GLUON DISTRIBUTION 



39 



A Momentum sum rule and gluon distribution 

In this appendix, we shall give the expression of the constant in the gluon distribu- 
tion, constrained by the momentum sum rule. Recall that we have, at = Qq, 

2 

xuv{x) = —x'^{l + jux){l-xy^, 



xdvix) = -^x''(l + 7dx)(l 

xqi{x) = [Alog^{l/x)+Blog{l/x) + C + Dix'^]{l-x] 
xg{x) = [Ag \og\l/x) + Bg log(l/x) + Cg] (1 - x)'+\ 



(5) 



where Ng is given by equation (2). Wc shall use momentum conservation to constrain 
the constant term Cg in the gluon distribution. Let us first introduce the special 
functions that we need. The Euler Gamma function is defined by 



r{x) = [ 

Jo 



oo 

dte-^e-K 





We can then introduce the Beta function B(x,y), the digamma function ^'(x) and 
the polygamma function *("*)(a;) related to the gamma functions by the following 
formulae 



T{x + y) ' 
d.T{x) 



r(x) ' 
^M(a;) = d^^{x). 

With these definitions, the momenta carried by the distributions (5) are given by 
the following expressions: 

2r? A ^ v + l \ A ^ V 



bu + ri + l\ '"6„ + r/ + 2y V ' 6« + r/ + l^ 

Pdv = 1. , , -I ( 1 + 7d I. , „ , n M 1 + 7d 



bd + r] + l\ '"6d + ry + 2y V 6^ + r? + 1 

2 



Pg, = ^ (^a|[7^; + ^'(6 + 2)]'-*«(6 + 2) + ^|+5[7£ + *(6 + 2)] + C^ 
+ DiB{b + l,rj + l), 

P9 = ^^(Ag!^[^E + ^{bg+2)f-^^'\bg + 2) + ^^+Bg[jE + ^ibg + 2)]+Cg 

Prom the proton, momentum conservation gives 

Pg+Puv +Pdv + 2iPu+Pd+Ps) = 1, 
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and we finally obtain 

+ {bg + '^■) - Puv - dPv - ^ {Pu + Pd + Ps)] ■ 
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